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We review models which generate a large non-Gaussianity of the local form. We first briefiy 
consider three models which generate the non-Gaussianity either at or after the end of inflation; the 
curvaton scenario, modulated (p)reheating and an inhomogeneous end of inflation. We then focus 
on ways of generating the non-Gaussianity during inflation. We derive general conditions which 
a product or sum separable potential must satisfy in order to generate a large local bispectrum 
during slow-roll inflation. As an application we consider two-field hybrid infiation. We then derive 
a formalism not based on slow roll which can be applied to models in which the slow-roll parame- 
ters become large before infiation ends. An exactly soluble two-field model is given in which this 
happens. Finally we also consider further non-Gaussian observables; a scale dependence of /nl and 
the trispectrum. 



I. INTRODUCTION 

There are many models of the universe which can predict a large non-Gaussianity. However the predicted amplitude 
and the shape of the non-Gaussianity are different among different classes of models. One category is those which 
generate the non-Gaussianity due to non-trivial classical dynamics on superliorizon scales. These models predict the 
shape of the bispectrum to be of the so called 'local type', which can be expressed as an expansion of the Bardeen 
potential P 

$(X) = $i(x) + /NL($i(x) - (<I>i(x))), (1) 

where <t> is the curvature perturbation on a Newtonian slice and $l is its linear and Gaussian part. (<f>|^(x)) denotes the 
ensemble average in a statistically homogeneous distribution. The current limit on the local type of the non-linearity 
parameter /nl from seven years of WMAP data [H is — 10 < /nl < 74 at the 95% confidence level. Constraints are 
expected to improve rapidly and significantly, first with Planck data and later using large scale structure data, see the 
recent reviews ■ The Bardeen potential is related to the primordial curvature perturbation of on large scales 
and in the matter dominated era by $ = (3/5)C. 

The curvature perturbation at horizon exit is determined by the classical perturbations of the scalar fields, 6(j)i{'x.). 
The subsequent evolution of ( can be conveniently described by the SN formalism (6l-[l0j . The curvature perturbation 
is given by up to quadratic terms [l^ 
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(2) 



where Af(x, t) is the c-folding number evaluated in an unperturbed Universe, from the epoch of horizon exit to 
later epoch of uniform energy density hypersurface (for an extension to include gradient terms see [ll|). The power 
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spectrum and the bispcctrum Bq are defined by 

(CkXk.) ^ (2^)3<53(ki+k2)^n(fci)> (3) 

(CkiCk.Ck3> = (2^)'<53(ki+k2+k3)Bc(fci,fc2,fc3). (4) 

From this we can define the observable quantities, the spectral index, the tensor-to-scalar ratio and the non-linearity 
parameter: 

_ aiog-Pc 
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where Pt = SVt^/Mp = 8H'^ / {An'^ Mp) is the power spectrum of the tensor metric fluctuations. It is well known 
that single-field inflation does not lead to a detectably large non-Gaussianity, in fact /nl is suppressed by slow- 
roll parameters [l^ . Observably large non-Gaussianity can be obtained by breaking the slow-roll conditions during 



inflation 13 1 using extended kinetic terms [ij] , see also the reviews [15|, |16[ , or going beyond models of single-field 
inflation [l3li3,[iill 

It is natural to consider multiple scalar fleld since they are ubiquitous in many beyond the standard model of particle 
physics, such as supersymmetry and string theory. These scalar fields generate non-adiabatic perturbations during 
inflation and change the evolution of the curvature perturbation after horizon exit. The residual isocurvature pertur- 
bation may be present in the primordial density fluctuation and can be correlated with the curvature perturbation 
or may be responsible for an observably large non-Gaussianity in the cosmic microwave background and large scale 
structure, for observational limits on isocurvature perturbations see (ll [23l. [23|. In this review we will only consider 
models with adiabatic primordial perturbations, in which the isocurvature perturbation present during inflation is 
converted into an adiabatic perturbation. We also neglect the secondary non-Gaussianities generated at later times, 
for example see [1^, . 

There are popular multi-field models (27j which may generate observably large non-Gaussianity. These include the 
curvaton scenario, modulated (p)rcheating and an inhomogencous end of inflation, sec Sec.|TTl In these scenarios, large 
non-Gaussianity is generated either by the means of ending inflation, or after inflation. It was shown recently that it 
also possible to generate large non-Gaussianity during the evolution of slow-roll multi-field inflation, see Sec. IIIII 

All of these models generate the large non-Gaussianity after horizon exit, such as after reheating, at the end of 
inflation, at the phase transition or during inflation after horizon exit and involves the perturbation of the non- 
adiabatic mode. Therefore the non-Gaussianity of these models is of the local type which is distinguishable from 
other shapes of non-Gaussianity (for a list of possibilities see e.g. [28l.[29|). in which the non-Gaussianity is generated 
intrinsically from the quantum fluctuations, or during horizon exit. 

In Sec. [TT]we summarise the aforementioned three models, which are popular methods of generating a large non- 
Gaussianity. In Sec. mil we review the possibility of generating a large non-Gaussianity from multi-field slow-roll 
inflation and in Sec. IIVI we consider hybrid inflation with two inflaton fields as an application. Then in Sec. |V] we 
discuss multi-field models of inflation without assuming the slow-roll conditions and present an exact solution. Non- 
Gaussian observables beyond /nl, such as its scale dependence and the trispectrum are introduced in Sec. lVIl Finally 
we conclude in Sec. I VIII 



II. SUMMARY OF POPULAR MODELS GENERATING A LARGE LOCAL NON-GAUSSIANITY 



Multi-field infiationary models [27[ can lead to an observable signature of non-Gaussianity. These include models 
in which the large non-Gaussianity is generated either by the means of ending infiation, or after inflation. We review 
the three most popular of these models in this section, for a discussion on how these models are related see js^]. In 
the next section we review the possibility of generating large non-Gaussianity during multi-flcld slow-roll inflation. 
Later in Sec. |V] we discuss multi-field models of inflation without assuming the slow- roll conditions. This does not 



exhaust all of thepossible ways of generating a large local non-Gaussianity from inflation, see also the reviews I3ll-|34 
and the papers [s^, . Non-Guassianity of the local form can also be generated in the ekpyrotic scenario |37l. l38l|. 
although in the simplest case /nl is large and negative [3^, which is observationally ruled out. 
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A. Curvaton scenario 



111 this scenario there is a Hght weakly-interacting inhoniogeneous scalar field, the curvaton in addition to the 
inflaton field, 4>. The additional scalar field is completely subdominant during inflation, however it can dominate the 
energy density of the Universe later since the time-averaged equation of state becomes that of the pressureless matter 
(P=0), whose energy density decreases more slowly than that of the radiation which the inflaton field has decayed 
into. After decay the inhomogeneity of the scalar field leads to the density perturbation of radiation [13, El- This 
is the basis of the curvaton scenario p^ - l43 | which liberates the inflaton field from being required to generate the 
observed spectrum of perturbations (although there is instead an upper bound on their allowed magnitude) . 

The significant non-Gaussianity of the local type can be generated in the curvaton model [l8| . A study of non- 
Gaussianity in the curvaton scenario using second-order cosmological perturbation theory was done by Bartolo et 
al [l9l | with a sudden decay approximation and also by Malik et al j2G | using a fully numeric approach, which ac- 
counted for both the sudden decay and non-sudden decay approximations. Lyth and Rodriguez O used the 5N 
formalism including the quadratic terms (second order terms in Eq. ([2])) to calculate the non-Gaussianity. The 
non-linear generalisation of the curvature perturbation and a numerical study was done to give a full probability 



distribution function in |45| . This simple curvaton model was extended into the mixed curvaton- inflaton scenario in 
which fluctuations of both the inflaton and a curvaton contribute to the primordial density perturbation studying 
the isocurvature perturbation in j46| . The multiple curvaton fields which contribute to the primordial density per- 
turbation was studied in (47l - l49j . Observational constraints on curvaton models with isocurvature perturbation and 
non-Gaussianity can be found in j50l - l52| . 

The generation of large non-Gaussianity in the curvaton scenario can be easily understood using the result of the 
non-linear SN formalism [IBl, [4^. Assuming the sudden decay of the curvaton field on the H = T, i.e., when the local 
Hubble rate equals the decay rate for the curvaton, uniform energy-density hypersurface leads to a non-linear relation 
between the local curvaton density and the radiation density before and after the decay, 

^2^g4(Cw-C)+f|^g3(c,-C)^l^ (8) 

where fir + ~ 1- Here C is the primordial curvature perturbation which remains constant on large scales in the 
radiation-dominated era after the curvaton decays and is the curvaton perturbation. Expanding this equation 
order by order yields up to second order 145} 

C = (1 - r^)Qnf + r^Cx + !:^ii^-I^i^±l2^(C^ - C„.f? + higher - order terms, (9) 
where r-^ = SQ-^/ (4 — fi^) at the time of the curvaton decay. The curvature perturbation of curvaton field [45| 

i/^y ao) 

3 xo 3 V Xo / 9 V xo / 

is almost constant before the curvaton decays from the start of curvaton oscillation. Note that the intrinsic non- 
Gaussianity in the curvaton field during oscillation is fNL,x ~ ~f i which comes from the non- linear relation between 
^iX and C^. After the curvaton decays into radiation, the curvature perturbation of radiation has the following 



aton decays 



non-Gaussianity parameters: 

/nl 



(l + f)2 
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where is evaluated when the curvaton decay based on the sudden decay approximation but gives good agreement 
with a full numerical study [10, HH] , and f is defined as 

" l(l-^x)C»/P- ^ ' 

In the limit of the pure curvaton scenario, Qnf — > 0, it recovers the usual normalization for /nl, i-c. ?^^/(l + '~)^ ^ 1- 
Therefore the large /nl is obtained for small r^. This large /nl, in spite of the order of unity non-Gaussianity in (^^, 
has the origin in the non-linear relation between (or Six) and C in Eq. 

In fact Eq. (|8]), and thus Eq. ([9]), can be used "at any time" before the curvaton decay with evaluated at that 
time to find the total curvature perturbation at that moment, as well as at the epoch of curvaton decay. This enables 
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us to understand the evolution of the curvature perturbation ( on the uniform energy density hypersurface before the 
curvaton decay. Initiahy when the primordial radiation was generated, e.g. from the inflaton decay, r^^ is negligibly 
small so we cannot ignore the perturbations of the inflaton field. At this early time, C — Qnf and /nl = since 
f (X which makes /nl negligible in Eq. (jlip . As time goes on increases while C,inf and Cx almost constant 
since there is no energy flow between two fluids satisfying adiabatic condition except around the decay time, thus the 
curvature perturbation evolves gradually. If the curvaton survives long enough, reaches a value with which the r^C^ 
dominates the first term, Cm/ , and the curvaton perturbation is then responsible for the primordial inhomogeneity of 
the Universe. The large non-Gaussianity is possible only when at the decay time is small, i.e. the energy density 
of curvaton is subdominant compared to that of radiation. This is because when = 1, the curvaton energy density 
dominates over the radiation component before decay, then ^ is effectively the curvature perturbation of a single field 
(curvaton), and therefore non-Gaussianity /nl = Inl,x = ~f • We will also see this phenomena in the multi- field 
slow-roll inflation later, in Sec. IIIII 



B. Modulated (p)reheating 

After inflation the energy density in the inflaton field must be transferred into radiation. In the simplest case 
of adiabatic perturbations this process does not affect the primordial curvature perturbation on scales which are 
observable today, because these scale were much larger than the horizon at the time of reheating. However in a 
spirit similar to the curvaton model, there may be a subdominant light scalar held present during inflation which 
modulates the efficiency of reheating. This makes the efficiency of reheating a spatially dependent process. The 
quasi-scale invariant perturbations in this field, which during inflation are an isocurvature perturbation, may be 
converted into the primordial curvature perturbation during this process. For a review of reheating after inflation see 
for example [H, [Hj] . 

As a simple illustration of this, in the "old" models of perturbative reheating the decay rate of the inflaton is given 
by r ~ A^m, where A is the strength of the coupling between the inflaton and reheating fields and m is the mass of 
the inflaton. Provided that T is much less than the Hubble rate at the end of inflation, so that reheating takes place 
slowly, the reheating temperature is given by 



Troh VrMp - X^/mMp . (13) 

If the coupling strength is a function of the local value of a second scalar held x then this will give rise to a perturbation 
in the decay rate of the inflaton field and thus in the reheating temperature which is responsible for the density 
perturbation after reheating. Provided that the x field is perturbed, the curvature perturbation can be written as 

C = , (14) 

where a is a parameter depending on the ratio of F to the Hubble expansion at the end of infiation, whose value 
increases with decreasing T/H and leads to 1 /6 in the limit F <^ Hi^f [U, [s^ . Hence perturbations of the light scalar 
field X can be imprinted into the radiation temperature through the reheating process [s^ [H^ . In a similar way the 
mass of the decaying particle may be modulated [s^l. The above equation can be extended to second order, which 
is required to calculate /nl- In the simple quadratic case that F = Fq -|- ^i{x/x)^ = Fq -|- Fi(l -|- Sx/x)'^ the result 

is [2il,[5i-i3 

5 Fo + Fi 

/NL = -r— f; ■ (15) 

Hence we can see that the non-Gaussianity can easily be larger than unity in this model by two sources. One is due to 
the small value of a which occurs when the decay rate of the inflation is not much smaller than the Hubble constant 
at the end of inflation. The other is when Fi is small in which x field only controls part of the decay channels of 
inflaton fleld. . 

It is now generally considered that a period of rapid and highly non-linear preheating preceded the perturbative 
reheating phase In this case one can instead consider modulated preheating, see e.g. [6ll - l63j . A particularly 

rapid form of preheating, known as instant preheating [63 | has attracted more attention in this area, as it is possible 
to find analytic estimates j65H69| . It is possible to generate the primordial curvature perturbation in these models by 
having either a weak (65j or a strong pl\ symmetry breaking of the two-field potential during inflation which affects 
the efficiency of instant preheating, or by having two inflaton fields which are coupled with different strengths to the 
preheat field However all of the above models of modulated instant preheating are only allowed, even at the 
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linear level in perturbations, for a small corner of parameter space and even then only for a special choice of initial 
conditions [69| . 

C. Inhomogeneous end of inflation 

An inhomogeneous phase transition can occur also at the end of inflation between inflationary and non-inflationary 
expansion due to the condition of how inflation ends [70| . In the case of multi-component inflation there is a family 
of inflationary trajectories. If the family of trajectories is not perpendicular to the surface on which inflation ends then 
there is a possibility that inflation ends on a slice which is not of uniform density and the generation of the curvature 
perturbation at the end of inflation is expected (22| . ^ In single field inflation the energy density is determined solely 
by the single field (j) and the end of inflation is parameterised by the unique value 0e- If there is another field a, then 
(peicr) will depend on position through the perturbation 5(t(x). As a result, the change of e-folding number from a 
spacetime slice of uniform energy density just before the end of inflation to a spacetime slice of uniform density just 
after the end of inflation will have a perturbation SN^ = Ce- This is the curvature perturbation generated by the 
end of inflation and can dominate over the contribution from inflation [2^ . 

Using the perturbation of 



(16) 



the curvature perturbation generated at the end of inflation is [22 



{5af 



(17) 



When (^e dominates the perturbation generated during the inflation, the power spectrum is given by 



a' 2 



2ee V 27r 



(18) 



The non-Gaussianity can also be generated at the end of inflation since the perturbation at the end of inflation is not 
required to be Gaussian. Thus we expect that the observable large non-Gaussianity can be generated in this process. 
From Eq. (|17p the non-Gaussianity is [l^ 



NL 



(19) 



This was applied to the simple example of extended hybrid inflation using sudden end approximation [7J, |75| . 

Recently Sasaki et al calculated analytically the curvature perturbation and the non-Gaussianity in the multi-brid 
inflation models [zl, [T^I from horizon exit to the end of inflation considering the general couplings to the waterfall 
held. In general we cannot just ignore the curvature perturbation generated during inflation. Furthermore that is also 
affected by the end of inflation. A deeper understanding of non-Gaussianity both from the evolution during multi-fleld 
hybrid inflation and effects from the end of inflation is considered by Byrnes et al [tsI [79j which is summarised in the 
next section. 



III. NON-GAUSSIANITY DURING SLOW-ROLL INFLATION 



It is well known that the single fleld slow-roll inflation with canonical kinetic terms generates a non-linear parameter 
of the order of the spectral tilt of the spectrum and thus too small to be observed [l2| • On the other hand, in multi-fleld 



^ This was also discussed in using the hnear relation between two fields at the end of inflation and the non-Gaussianity is transferred 

from that of the light field which is generated from its non-linear self-coupling. 
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inflationary models the non-adiabatic perturbation may change this conclusion. Furthermore even within slow-roll 
inflation an observably large non-Gaussianitycan be generated if certain conditions are satisfied. 

Rigopoulos et al introduced a formalism [sO, HH to deal with non-linearity based on a system of fully non-linear 
equations for long wavelengths and applied it to a two-field inflation model. In the case of two-field quadratic inflation 
with two massive fields, an analytic calculation is possible with slow-roll conditions and this shows that the nonlinear 
parameter is small (sH, . The numerical calculation confirms the analytical prediction that it is possible to generate 
a narrow spike of large non-Gaussianity while the inflaton trajectory turns a corner, but the non- Gaussianity decays 
quickly after the corner (ssl. Is^. This spike is due to the temporary jump in the slow-roll parameters [ssl ISGj. 

Recently an analysis based on the general analytic formula was done by Byrnes et al [tI, [t^ and they showed 
that it is possible to generate observable large non-Gaussianity even during slow-roll multi-field inflation models and 
explicitly wrote the conditions as well as concrete examples^. In this section we summarise the general conditions 
for two-field slow-roll inflation to generate a large non-Gaussianity and briefly give a speciflc example. In the next 
section we consider the multi-fleld hybrid inflation model which has been studied in depth and we include the effects 
from the end of inflation. 



A. General formulas 



In the case of a potential which is either sum or product separable, it is possible to compute the curvature perturba- 
tion using the (5A^-formalism and slow-roll conditions, and an analytic formula for the non-linearity parameter can be 
obtained. The general formula for the non-linearity parameter /nl was calculated in the case of a separable potential 
by sum 18 31 and this was later extended to the non-canonical kinetic terms as well as to the separable potential by 
product |88j . This has been further generalized to the arbitrary number of fields [13 and to the trispectrum [oO] ■ 

For a product separable potential, W{ip^x) = U{ip)V{x)j can find the number of e-foldings analytically in the 
slow- roll limit [qi} . 

Then the power spectrum and the non-linear parameters can be calculated straightforwardly using the derivatives 
of e-folding number [ssl Issl . [9lj. In the case of two inflatons with canonical kinetic terms, the power spectrum and 
spectral index are found to be [s^ 

W. (u' if\ , , 



— 1 = 



-2e* -4- 



2e' 



''XX 

2et 



(22) 



r = 16(^ + ^) , (23) 



where the slow-roll parameters are 



^^^^^-V.,eos^,, .^M1(Y^] .,3in2., (24) 



where 9 is the angle between the adiabatic perturbation and one of the flelds and 



_ n .7-2 ^vv A j2 Uy^x — m'^Ym. (2F>) 



^ There is a claim [stII that the sizeable value of /nl can arise from a loop corrections. 
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Note that the superscript "*" and "e" denotes that the values are evaluated at horizon exit and at end of inflation 
respectively. In the above u and v are evaluated at the end of inflation as 



U = — = COS I 



v=^= sin^ 61"=. 



(26) 



The non-linearity parameter /nl becomes [8S 



where 



/nl — ^ 





e* 



At 



2e; 



e e 



(27) 



(28) 
(29) 



Similar formulas are found in the case of a sum separable potential j83l . ISS 



B. Conditions for generating a large /nl 

Using an analytic formula for the nonlinear parameter /nl in the case of a sum or product separable potential, 
it is possible to generate significant non-Gaussianity even during slow-roll infiation with Gaussian perturbations at 
Hubble exit and the general conditions for it can be written explicitly [t^ . 

For the separable potential by product large nonlinear parameter |/nl| ^ 1 is possible when 



sm 



2 n* 



< 



sm 



4 ne 



1 



- 1 



/sin^ e^Gp 

With this condition, the large /nl is given by 



= 5 \-^xx 



2t] 



XX I 



(30) 
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NL 
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2< 



(31) 



For large /nl the substantial increase in sin 6 is necessary between horizon exit and the end of infiation, which 
corresponds to the curve of the trajectory in the field space during the slow-roll infiation. However the background 
trajectory must be almost entirely in the one-field direction. 

One specific toy model is the quadratic times exponential potential [tsI ] where the potential is given by 



(32) 



With expansion of exponential, this potential can be understood as the quadratic chaotic infiation model which has 
a coupling to another light scalar field. For this example large /nl is generated when the (p at horizon exit is small 
ip* ~ IQ-^Mp for A = 0.04 and x* = 16Mp. 

For a sum separable potential, we can find similar conditions for generating a large /nl [ill ■ 



^ There is a symmetric region of large non-Gaussianity under the exchange of the two fields. 
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IV. HYBRID INFLATION WITH TWO INFLATON FIELDS 



We consider a model of two field hybrid inflation, whose potential is given by 

= Woexp j exp (^2'?xx^ + -rj,,j^ + j , (33) 

which is vacuum dominated, i.e. which satisfies |?/,^,^<y5^| <Si Mp and j^xX^I ^ Mp. It is in this regime that the two 
ways of writing the potential given above are equivalent (at leading order in slow roll), so we can treat the potential 
as being both sum and product separable. We assume that inflation ends abruptly by a waterfall field which is heavy 
during inflation and hence doesn't affect the dynamics during inflation. First we calculate observables during slow-roll 
inflation. We will consider the full potential including the waterfall field in Sec. IIVBI which incorporate the effects 
from the end of inflation considered in Sec. Ill CI We will see that this can lead to a change in observables on the 
surface where the waterfall field is destabilised. 

In the vacuum dominated regime the slow-roll solutions are 

ip{N) = v'^e-"-^, = X^e-"-^, (34) 

where denotes the value at the horizon exit. Throughout this section whenever we write a quantity without making 
it explicit at which time it should be evaluated, we mean the equation to be valid at any time N e-foldings after 
Hubble exit and while slow roll is valid. Generally we will be interested in quantities at the end of infiation, in which 
case wc take = 60. 

The slow-roll parameters are 

I 2 I 2 X^ 

= o^mmTtt, '^x ^ o'yYvTTT' e^t^p + t^. (35) 



We note that the dominant slow-roll parameters r]^^ and r]^-^ are constants during infiation in the vacuum dominated 
regime and that they are much larger than the slow-roll parameters tip and throughout inflation. 
From the previous section, large non-Gaussianity can be realised in either of two regions 

cos20=^/^^^^«f, or sin20=— ^~^^«f. (36) 

V + X + Ex V + X + ex 

Since the two regions are symmetrical (TSj (before specifying the values of ri^p^p and fy^x)' "-"^ ^^^^ section we 

will focus on the second region. In this region where ^ e^, |/nl| > 1 is fulfilled by the condition. 




in other words. 



This implies three inequalities on B: 



sin^r<sin^^h/^-I , (37) 



'yxxr'e-^^"--"-)^ < sin^ ^ < hxxl- (38) 



2.* ^ 1 (^\^ (^\\ ,2 .2. 5 , , sin^fl 24 1 



Note that in this region sinfl ~ %xX/(^v"^'^)' fro™ Eq. p4p we require N[riipip — '7xx) > 1 so that siv? grows 
significantly during infiation. 
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A. Simplified formula for the observables when /nl is large 

We can substantially simplify all of the above formula in the case where /nl is large. We define the quantity 

2 / n AT \ 2 



r ; 



/dN\ ,(dN\ ^ fxg2(,,,-,,,)7v (40) 



In the region we arc considering where /nl is large, this is approximately given by the initial and final angles of the 
background trajectory with different exponents 



sin'' 6 , , 

-T^- (41) 
sm {)* 



In the case of large non-Gaussianity it follows that 



Vr { 1 + ^2^e2('^^^^-''xx)JV^ = ^ /£l V (42) 



„ 1 ^ 2^^^^^tl!Z2« , (43) 
I + r 



sin*^ 6" 5 



The first condition in Eq. (|39p implies that 



X*<</3*. (46) 



We therefore require a very small value of x* in order to have a large non-Gaussianity. While this may be a fine 
tuning, we comment that requiring a large ratio of the initial field values is common in models where the isocurvature 
perturbation generates a significant part of the primordial adiabatic perturbation, see the discussion in the conclusion 
of [69j . If the infiaton perturbations are neglected a priori, such as in the pure (i.e. not mixed) curvaton scenario, 
this remains implicitly true. For an explicit discussion of this, in the simple case that both the infiaton and curvaton 
fields have a quadratic potential see [92l ■ 

The sign of /nl is determined by the sign of 77^^ . The amplitude of /nl depends exponentially on the difference 
of the slow- roll parameters, rj^p^p — rj-^-^, which we require to be positive to be in the branch of large non-Gaussianity 
where sin^ <C 1, while the spectral index depends on a weighted sum of the slow- roll parameters, so it is possible to 
have a large non-Gaussianity and a scale invariant spectrum. However it is not possible to have a large and positive 
/nl and a red spectrum of perturbations. We will see in Sec. lIVB^ that by including the effect of the waterfall field 
this conclusion may change, depending on the values of the coupling constants between the two infiaton fields and 
the waterfall field. 

In Table U we give some explicit examples of values of ?7ipip, ??xx' V'* x* which lead to a large non-Gaussianity. 
Using Eq. ([S]) we also calculate the spectral index. The first example in the Table U shows that it is possible to have 
|/nl| — 100 and a scale invariant spectrum. We also see that it is possible to generate a large non-Gaussianity during 
slow roll with 77^^ and Ty^x both positive or both negative, or when one is positive and the other negative corresponding 
to a saddle point. The results for this model were verified using a novel calculational method of momentum transport 
equations in (93| . 



B. Effect of the waterfall field and further evolution after inflation 



In this section we include the effects of the waterfall field p which is required to end hybrid infiation. Infiation ends 
when the waterfall field is destabilised, i.e. when its effective mass becomes negative. During infiation the waterfall 
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-11.1 


357 


-2.6 


-0.02 


0.002 


-0.01 


-0.09 


1 


3 X 10-6 


0.16 


-132 


1.8 X 10^ 


-44 


-0.04 


0.0007 


0.06 


-0.01 


1 


4.3 X 10^^ 


0.1 


-3 


148 


-0.2 


0.11 


0.026 


0.01 


-0.06 


1 


7.5 X 10-6 


0.04 


-8 


2.5x10^ 


-2 


0.01 


0.0008 



TABLE I: Table showing some initial conditions for the hybrid inflation model that lead to large levels of non-Gaussianity. 
The table shows the parameter f, the bispectrum and the trispectrum non-linearity parameters, the spectral index and tensor- 
to-scalar ratio. They are evaluated when the number of e-foldings from the end of inflation is Nk — 60. The trispectrum 
parameters tnl and qnl are defined and calculated in Sec. IVI Al 



field is heavy and it is trapped with a vacuum expectation value of zero, so we can neglect it during inflation. The 
end of inflation occurs vifhen the effective mass of the vifaterfall field is zero, which occurs on a hypersurface defined 
in general by [zl, [t^] , 

a2 = G(^,x) = .gV+.92V, (47) 
which is realised by the potential W{(p, x), defined by Eq. ([55)1 . where Wq is given by 

Wo = Ig{^,x)p' + ^(^p'-^)' ■ (48) 

Here gi (32) is the coupling between the ip (x) field and p is the waterfall field with self coupling A. In general the 
hypersurface defined by this end condition is not a surface of uniform energy density. 

As discussed earlier, this is an example of a model with an inhomogencous end of inflation, i.e. where inflation ends 
at slightly different times in different places. It has also been shown for the hybrid potential we are considering that 
this can be used to generate a large amount of non-Gaussianity, for certain parameters values and fine tuning of the 
parameters [tsI [77| . However these papers concern the large non-Gaussianity generated at the end of inflation rather 
than during slow-roll inflation, by having a very large ratio of couplings (71/52 ^ 1- Here we consider the case where 
gi and 52 have the same order of magnitude with g\/g2 = Vw/Vxx sec. lIVBTl and with g^ = 52 in sec. IIVIT2] 



In this case we have chosen the coupling constants (which can satisfy gi <Q and/or < 0) such that the surface 
where the waterfall field is destabilised corresponds to a surface of uniform energy density. This is because the end 
condition can be rewritten as cr^ = 2g2/ri^x(h''lvv^'^ + I^xxX^) ^ In this case the value of all observable quantities 
such as the power spectrum and non-Gaussianity are the same as those we calculated previously which were valid at 
the final hypersurface of uniform energy density during inflation. 



2 2 
5i = 92 



In this case, the end of inflation given by the condition in Eq. (|47p docs not occur on a uniform energy density 
hypersurface [77|. We will show how the non-Gaussianity is modified by the condition at the end of inflation in this 
example. In general we expect there to be some modification to non-Gaussianity from the end of inflation, except in 
the special case we considered in IIVB II We have checked in [t^ that the correction from the extra expansion which 
occurs from the surface on which inflation ends up to a surface of uniform energy density is small and therefore one 
can still use the SN formalism in this case. 

In this case the power spectrum and /nl in the limit of large non-Gaussianity are [77j 



11 





f* X* 
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/NL 


Tl^ — 1 
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1 O.oXlU 


1 


1 OQ 
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u 


u.uuo 


04 -0 04 


1 1.5 X 10^"* 


5 


68 
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0.002 


0.08 0.01 


1 0.0018 


1 


4.59 


0.02 


0.0008 


0.02 -0.04 


1 0.00037 


1 


3.5 


0.02 


0.026 


-0.01 -0.09 


1 3 X 10-^ 


0.16 


-0.2 


-0.02 


0.0008 


0.06 -0.01 


1 4.3 X 10--* 


0.1 


38 


0.01 


0.006 



TABLE II: Same as Table |I] but with different end condition, = as used in sec. II V lT2l 



(50) 



16e* 1 



(51) 



/NL — T' 



2 '/w^ 



(52) 



We can see that the observables at the end of inflation are changed by the additional ratio of f]^y/?7^^ before f. 
Furthermore, /nl has an additional factor of rj^p^p/rj^^. 

In Table HI] we show the values of /nl, n/^ — l and r with end condition gl — g\ for the same parameter values that 
we used in Table H] where the inflation ends on the uniform energy density hypersurface. The first two examples show 
that if r]^^f, = —t\^x then the observables are unchanged for two different end conditions except that the sign of /nl 
is switched. The second example in the table shows that in this case it is possible to have a red spectral index and a 
positive value of /nl- For many values of the initial parameters when \'r\LpLp I f]xx\ ^ ^he magnitude of /nl decreases 
compared to Table U but the reverse can also happen, an example of this where |/nl| grows by more than an order 
of magnitude is shown in the final row of Tables U and [ill 



3. Further evolution after inflation 



So far we have assumed a quick transition to the radiation epoch at the end of inflation, thereby neglecting 
the dynamics of the waterfall field. However if we consider the role of the waterfall field, then after the waterfall 
field is destabilised there may be a further evolution of the primordial curvature perturbation, which will lead to a 
change of the observable parameters. This applies to any model with an inhomogeneous end of inflation since there 
are isocurvature perturbations still present after the waterfall field is destabilised and inflation has ended. Further 
evolution will depend on the details of reheating in a model dependent way. To the best of our knowledge this issue 
has not been considered in depth in any paper. If we assume an instantaneous transition to radiation domination (so 
a completely efficient and immediate decay of the waterfall and inflaton flelds) then there will be no further change 
to the observables as wc have argued in the previous section. However this is clearly an idealised case. 

In the special case where the waterfall field is also light during infiation Barnaby and Cline (sJ] have shown there is 
the possibility of generating a large non-Gaussianity during preheating for certain parameter values. This is possible 
even if there is only one inflaton field and the waterfall field present. However in this case inflation does not end 
abruptly when the waterfall field is destabilised so this is not the scenario we have considered in this paper. A recent 
claim that even a heavy waterfall field may lead to a scale invariant local non-Gaussianity has been withdrawn [951] . 



V. MULTIPLE-FIELD INFLATION WITHOUT SLOW ROLL 



In this section we demonstrate a new method that provides, for certain classes of models, analytical expressions for 
/jv L valid in regimes beyond a slow-roll approximation (96j . This method is based on the first-order Hamilton- Jacobi 
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formalism developed by Salopek and Bond [97|, which allows us to express inflationary observables in the multi- field 
case, without having to focus on a slow-roll regime (see for a similar application of this formalism to the single-field 
case). 

The main idea on which our arguments are based is to focus on infiationary multi-field trajectories in which the 
Hubble rate, and not the potential, is separable. This approach allows us to analytically study cases where the non- 
Gaussian parameter /nl becomes large, in regimes in which the slow-roll parameters are enhanced after Hubble exit, 
but before inflation ends. In particular, we provide analytical equations that express the non-linearity parameter /nl, 
in terms of quantities that generalise the slow-roll parameters during inflation. Using this method we demonstrate in 
the next subsection an exact solution of multi-field inflation which can give rise to a large non-Gaussianity due to a 
strong break down in slow roll shortly before inflation ends. 

The first-order Hamilton- Jacobi equations of motion are 



dH 



= -2- 



dH 



dH 
dx 



o dH 
dx 



We concentrate on solutions satisfying the following Ansatz for the Hubble parameter H((j), x)' 

H{cl>, x) = i?W(0) + i^^'Hx), 



(53) 
(54) 

(55) 



that is, we demand that it can be split as a sum ^ of two pieces, each one depending on a single field. Notice that 
this is different with respect to the requirement of separable potentials usually done in the literature and earlier in 
this article. 

It is convenient to deflnc the following quantities: 



r 



H 



(1)' 



H 



H 



(1) 



H 



H 



(2)- 



H 



H 



(2) 
XX 



H 



(56) 
(57) 



and 5 = 5'^ -\- 5^ . Although their deflnition resembles the corresponding one for the usual slow roll parameters e and 
77, they do not coincide with them when taking a slow roll limit ^ . We assume that the quantities 5 and 7 are much 
smaller than unity at Hubble exit t = t^, in order that we can use the 5N formalism. During inflation, by definition, 
we have to ensure that the quantity tn = —H/H^ = 25 <1. The quantities 7 can however become much larger than 
unity during inflation. 

Using similar techniques to those applied to a separable potential, one can calculate the spectral index and /nl, 
the details were given in [9^ and the results are 



riQ — 1 = —4(5, 



UH (1 - -^"-ff) + Wi/ (^1 - ^Wif) 



(58) 



/NL 



St 



1 1: 



7. 



2 ( 



Vh 



Ah 



(59) 



Hence wc follow an analysis related to ISal. In an analogous way, one could also consider a situation in which the Hubble parameter 
splits as a product of single-field pieces [8Sl | . The formula's in this section were extended to an arbitrary number of fields with a sum 
separable Hubble parameter in l9il . 

It is nevertheless simple to work out the relation with the slow roll parameters, defined by Eqns. H24|l and H25|l . In a slow-roll regime, 
one finds 

2 ' 2 ■ 
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Z,= [H(^'^St-Hi'^S^)/S,, (60) 

We reiterate that the formula for /nl is exact (apart from neglecting a slow-roll suppressed contribution due to 
the non-Gaussianity of the fields at Hubble exit) and is not based on a slow-roll expansion. Although 6^ and 7* are 
necessarily small, the quantity A is not suppressed by these parameters, and can assume large values during inflation, 
enhancing /nl ■ We note in agreement with [83| that if one of the fields has reached a minimum so = or % = 
at the end of inflation then Ah = and the non-Gaussianity at the end of inflation will be very small. Therefore 
any model of inflation with a separable potential or separable Hubble factor with a large non-Gaussianity present at 
the end of inflation must have both fields still evolving, and therefore the presence of isocurvature modes. It would 
therefore also be interesting to stud y th e evolution of the perturbations after inflation, to see if this leads to an 
important change in the observables jlOOj . 



where the new symbols arc denned by 



Exact solution 



Very few exact solutions in multiple field infiation are known jlOlj . Here we present an exact solution which can 
give rise to a very large non-Gaussianity at the end of i nfiati on. A potential with a similar form may be motivated in 
string theory, in the context of Kahler moduli inflation jl02l | . 

The potential we consider is [9^ 

W{4>, x) = C/o (^1 - Aie-""^ + ^se-'""^ - B^e-^^ + 626-^^^ + e-"^-'^^^^ , (62) 



where the parameters a and /3 satisfy 

^'^"7;-^. P'' = "t:^^- (63) 



2_3_6A2 2 _ 3 6^2 

2 Al ' ^ ~ 2 Bf 



The Hamilton- Jacobi equations have the following exact solution for H and the fields: 

H ^ Ho[l- ^ie-"^ - ^e-^^^ , (64) 

^ = i In [e"^- -A,a^H„t], x = ^ In [e^"" ~ B^ Ho t] , (65) 
a P 



where we have defined Hq = ■\/C/o/3, t is cosmic time and we have set = 0. The overall factor to the potential Uq can 
be freely chosen so that the amplitude of the scalar power spectrum matches the observed amplitude of perturbations 
in the CMB. The scale factor results: 

a{t) = ao (e"** - Ai Hq t) ^ (e^'^* - /^^ Bi Hq t) ^ e"°K (66) 

Notice that the solution becomes singular at late times, when the scale factor vanishes and the field values diverge. 
This singularity occurs well after inflation ends, and we will not need to discuss it in our analysis. From the relations 



= -aVS'f' , 7>^ = -/3V(5x, (67) 

we notice that in this example 7"^ and 7-^ can become much larger than unity at the end of inflation if |a| and/or 
are much greater than one. This corresponds to a break down in slow roll, although en remains smaller than unity 
during the inflationary era by definition. 

As a simple, concrete example of parameter choice which gives rise to a large non-Gaussianity, we choose a regime 
in which the parameters |q;| and \(3\ are both large, let us say larger than some quantity R ^ 1. We note from 
that we are therefore required to have A2 < and B2 < 0. In this case the potential does not have a minimum, but 
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inflation still ends through eu growing larger than unity, and we can trust our results in this regime. The potential 
would need modifications which apply after inflation in order for reheating to take place and these modifications may 
provide a minimum for the potential. We parameterise the values of the quantities 5f and 5^ at the end of inflation 
as 



5t = ^ , 5^ = such that — + - = 1 . (68) 
2m 2n m n 

where the final equality follows from e^f = 1. We assume that the quantities m and n are not too big. Namely, they 

satisfy the inequality m, n <^ R^. This implies that Hq ~ i7» ~ He- We choose \jt\ = |7?| = 4^, such that we 
can write 

Ntot = 60 ~ 10 (p + q) . (69) 

This fixes the initial values for the fields and x*, to the values 

e""^* = 20pa'^Ai, (70) 

e/^^' = 20 q 13'^ Bi. (71) 



The conditions ([68]) imply 



Hot, - 20p- J— , (72) 



Hote = 20?-^^. (73) 
Since m and n are much smaller than i?^, we expect p q. We hence have the simple relations 



*-mlta^' *-ANltP^' 2N 

From (l58l). we find 



nc-1 ^ - ^„0.04. (75) 



So we have a red spectral index in agreement with present day observations [103| and in j96| it was also shown that 
the tensor-to-scalar ratio for this model is negligible. 
When evaluated at the end of inflation, wc find 



5 



1 ^ 













(76) 



We stress that this formula provide only the dominant contributions to /nl- It is valid in the case |Q;|,|/3|>i?^l, 
and m,n < . 

As an explicit case, we take a — — 100,/3 = 20, m — ~ 6/5 and the initial conditions to satisfy ([Tn]) and ([7T|) 
with Ntot = 60. Then, independently of the values of Ai and Bi, wc find f^L — 58 from the simplified formula (|76p . 
and /at i ~ 53 from the full formula ([5^ . See figures [T] and [2J Notice that the results for f^L differ from each other 
by quantities of order 1/R = 1//3 as expected by the approximations we have made. At the end of inflation, the 
parameters are 7*^ = 29, 7^ = — 13 and this shows that the slow-roll approximation has been significantly broken, see 
footnote [51 From Fig. [5] one can see that the trajectory is straight for much of inflation but it turn near the end of 
inflation (during the last c-folding before inflation ends), and it is during this time that /nl grows larger than unity. 

Notice that our requirements of final values for the quantities 5f and 5^ (both much bigger than 1/i?) imposes 
fine-tuning constraints on the initial conditions, since at leading order in 1/i? the values for 7^ and 7? must coincide 
(see Eqs. ([7n|)-([7T]) and rccaU that p^q). 

We stress that the enhancement of non-Gaussianity occurs towards the end of inflation. What happens just after 
inflation is a model dependent issue; we cannot address this question within the approximations used in this concrete 
model. For our choice of parameters we notice that /at l is still increasing at the end of inflation, see figure [TJ but for 
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fNL 




FIG. 1: Plot showing /jvl as a function of eH towards the end of inflation, for the values of the parameters given in the text. 
Inflation ends when en = 1; for this example /jvz, ~ 59 at that time. 




FIG. 2: Left plot shows the trajectory considered for the parameters given after eq. (|76p superimposed on a contour plot of the 
potential. The square on the trajectory indicates a point along the trajectory one e-folding before inflation ends as (f) and x roll 
towards zero. This shows that the fields roll much more quickly during the final stage of inflation, and the trajectory curves 
near the end. The right plot shows the potential for the same parameter values. Notice that inflation ends on the plateau long 
before the potential becomes negative. 



other choices it may start to decrease before inflation ends, see [99|. We stress that in the regime where \a\, \/3\ ^ 1 the 
slow-roll parameters will necessarily become much greater than unity by the end of inflation, which may correspond 
to 17*** I ^ 1: this is one reason by which f^L can become large. Our formalism is at least in principle suitable to 
study the evolution of non-Gaussianity after inflation, when the parameter en becomes larger than unity (although 
we restate that the specific potential being considered here anyway needs to be modified after inflation ends) . 

VI. HIGHER ORDER NON-GAUSSIAN OBSERVABLES 

Non-Gaussianity is most popularly parametcrised in terms of one non-linearity parameter, /nl- There are several 
different /nl's used in the literature which parameterise different shape dependences of the bispectrum [2^, the most 
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popular being the local model discussed in this article and the so-called equilateral model which can arise from models 
of inflation with non-canonical kinetic terms (a popular example is DBI inflation [31 )■ But can we learn more than one 
number from an observation of non-Gaussianity? And if /nl**' is detected how can we hope to distinguish between the 
models which can generate this? Fortunately observations of non-Gaussianity have the possibility of teaching us much 
more than the value of one parameter. Two examples of further non-Gaussian observables are a scale dependence of 
/nl and the trispectrum. 



A. Trispectrum 

The four point function of the primordial curvature perturbation is defined by 

(Cki Ck. Ck3Ck4>e = rc(ki,k2,k3,k4)(2^)353(ki + ^2 + ks + U) , (77) 

which using the SN formalism, and assuming that the fields have a Gaussian distribution at Hubble exit, is given by 

Tc(ki,k2,k3,k4) =TNL [Pciki3)Pc{k3)Pc{ki) + (11 perms)] 
54 

+ ^9NL [Pc{k2)Pc{k3)P^{ki) + (3 perms)] , (78) 



where ^13 = |ki + k3|, and the trispectrum non-linearity parameters can be calculated using [82l. l90l llO 

25 NabcN'^N^N^ ^^^^ 
- 54 (NoN^r ■ ^ ' 

Hence we see that the trispectrum depends on two non-linearity parameters (as opposed to one, /nl, for the bis- 
pectrum), and they may be observationally distinguishable since they are prefactors of terms with different shape 
dependences (jTS]) . The current observational bound on the local type of the bispectrum from seven years of WMAP 
data is — 10 < /nl < 74 at the 2a level Recently there has been a first constraint on both t^vl and gNL: which cam e 
from using WMAP5 data, with the bounds -7.4 < 5jvl/10^ < 8.2 and -0.6 < t^vl/IO"* < 3.3 at 95% confidence p^ . 
There have been two other observational constraints on the trispectrum through gNL, setti ng tnl — 0. The bou nds, 
which in both cases are roughly \gNL\ ^ 10^ — 10^, come from large scale structure |l06| . and the CMB [l07| . If 
there is no detection of non-Gaussianity it is expected that with Planck data the bounds will be reduce d to abou t 
|/nl| ^ 10, TjvL ^ 10"^ and g^L ^ 10^ at the 2a level and future observations may become even tighter |l08l . 
Note that the expected future constraint on gjvL is about two orders of magnitude weaker than that on t^vl 
We will see that it is also possible for some models and parameter ranges that the first observational signature will 
come through the trispectrum rather than the bispectrum. 



m 

IOC 



1. Trispectrum in two-field hybrid inflation 
In the regime where \gNL\ and tjvl are greater than unity they are given by 



_ 10 ^(^yy"2??xx)-??xx 4- /SIN 

gNL - ^ ——: /nl, [oi-) 

3 1 + r 



tnl 



{l + f) 



f \3 'XX 



2 



" e^^^"--"-) = ( I/nl ) ■ (82) 



We see that g^L is subdominant to /nl and hence won't provide a competitive observational signature. It follows 
from that t^l > (6/nl/5)^, so tml may be large and provide an extra observable parameter for this model. This 
inequality between t^l and /nl is true in general [59[ , and equality is reached whenever a single field direction during 
inflation generates the primordial curvature perturbation. However it is usually assumed that tjvi, ~ /^l since both 
arise from second derivatives in the 5N formalism. In fact for our model it is possible to have a small /nl (and hence 
also a small gNL) but a large and potentially observable tnl- For this we require that f ^ 1, although in practice 
if we make it too small it may no longer be possible to satis fy a constraint on a minimum possible field velocity 
consistent with a classical slow-roll trajectory, as discussed in |79| (see also |llOj and a discussion which reaches a 
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different conclusion is given in [llll ] ) . In the final example in Table |T] we give an explicit example of parameter values 
which give rise to an /nl which is probably too small to be detected with Planck but with a very large trispectrum 
through tnl > 10'^ that should be detectable at a high significanc e. For another example with fNL,gNL ^ 0{1) but 
tnl ^ 1, see [gO]- In contrast it has been shown in several papers jl04l 1112141161 that in the curvaton scenario where 
the curvaton has a non-quadratic potential it is possible to realise \gNL\ 3> 1 while tjvl = (6/nl/5)^ is small with 
some tuning of parameters. This is also possible in the exact solution which we presented in the previous section and 
we consider this next. 



2. Trispectrum in the exact solution 

In this regime where \a\ ^ |/3| > i? and m,n > 1/R^ it is also possible to give compact expressions for the 
trispectrum (4-point function) non-linearity parameters, in terms of as 

/6 y 54 3n-m/6 V 

rNL-[^-fNLj , ^5iVL - -2^^ (^g/^Lj . (83) 

The complete result was given in [9^. Hence both of the trispectrum non-linearity parameters are generally large 
whenever /nl is. We note that if m ~ 1 then from (|68|) n 3> 1, so the trispectrum through g^^ will give the dominant 
signal of non-Gaussianity through a large, negative gNL- 



B. Scale dependence of /nl 

In its simplest form, the local form of /nl defined by Eq. ([1]) is a constant parameter, independent of both position 
and scale. However in realistic models /nl is likely to be mildly scale dependent. This may happen in two ways. 
In general single-field models (such as the curvaton scenario) where the scalar field which generates the primordial 
curvature perturbation has a non-quadratic potential the non-linearities this generates will give rise to scale dependence 
of /nl- Alternatively even in models where all of the fields have a quadratic potential; if the primordial curvature 
perturbation has contributions from more than one field, and the fields don't all have the same mass, then /nl will 
again have a scale dependence. This is because the correlation between the first and second order terms of C, which 
the bispcctrum depends on, will become scale dependent. This is indeed what happe ns in the two-field hybrid model 



presented earlier. For detail s about a scale dependence of local /nl see Byrnes et al jll7 |. 
We define a second observable parameter derived from the bispectrum 

which is analogous to the scale dependence of the power spectrum. There is a subtlety here, which is that in general 
/nl may depend on three independent parameters, fci,fc2 and fca. However it was shown in jll7j | that provided one 
takes the derivative while keeping the ratio of the three k vectors fixed then rifj^-^ is independent of the shape of 
the triangle described by the thre e k vectors, which makes this a well defined quantity. Observational prospects for 
this quantity were considered in [llSf . who showed that the Planck satellite is sensitive to nj^-^ ~ 0.1 assuming a 
fiducial value of /nl = 50. This observational sensitivity to 't./nl is about a factor of two larger than the current 
preferred value of the power spcctrums spectral index and may therefore provide an interesting extra constraint on 
non-Gaussian models. The scale dependence of an equil ateral for m of /nl has also been considered from both a 
theoretical and observational perspective, see for example [TTsHmt . 



In the example of two-field hybrid inflation /nl has a scale dependence both because of the exponential term in 
/nl, P5|). and because r will vary through the change of the initial value of sin^ 9* . We find 

dhTk ^ m^k = -^^""-^ 

Using this we find from (HSl) that 

nf..^^^^^--^^^. (86) 
dlogfc 1 + r ^ ^ 



In the case that we include the effect from the surface where the waterfall field is destabilised and gf = g| we find 
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from ([52]) that 



4 



(87) 



1 + (!Lp£.] 



r 



For both cases the spectral index of /nl satisfies 



Vxx) < 



NL 



<0, 



(88) 



for any value of f and hence /nl will be smaller on small scales. 

Because we require a relatively large value of 77^^ — ry^^ > 1/A^ for our model to generate a large non-Gaussianity 
it is quite possible for our model to generate a relatively significant scale dependence of /nl- However the amount 
also depends on f and when this is large then the x field is almost solely responsible for generating C, at both first 
and second order and Uf^-^ is suppressed. We can also see in agreement with the statement at the beginning of this 
subsection that in this case if the mass of the two fields are equal then the two fields have the same scale dependence 
and Uf^^ =Q. ^ ^ ^ ^ ^ ^ 

We note that this is in contrast to the large non-Gaussianity from an inhomogcneous end of inflation found in |77| . 
In the specific cases they considered to generate a large non-Gaussianity the non-Gaussianity was generated purely 
at the end of infiation and /nl is scale independent. In detail we see from Eqs. (4.4) and (4.24) in [tJ that their 
formulas for /nl does not depend on N or on any quantities evaluated at Hubble exit. This is also in contrast to the 
exact solution considered in section IV Al Our expression for /nl in eq. (j76p depends on the initial values through p 
and g, and this dependence drops out at leading order in a. Hence f^L is independent of the number of e-foldings 
and scale independent. 



We have reviewed various models which can generate a large local non-Gaussianity. A feature shared by all of 
these models is that they have more than one light scalar field present during inflation. This extra degree of freedom 
generates an isocurvature perturbation which is at least partially converted into the primordial curvature perturbation 
after horizon exit of the modes which are observable today. In the curvaton and modulated reheating scenarios this 
conversion occurs after the end of inflation, while in the inhomogcneous end of inflation scenario this conversion occurs 
on the non-uniform energy density hypersurface on which infiation ends. For these three scenarios the light field which 
generates the primordial curvature perturbation after or at the end of inflation can be treated as a test field which 
does not affect the inflationary dynamics. 

Our main focus has been on models in which a large non-Gaussianity is generated during inflation. This can occur 
even within slow-roll inflation for certain potentials and certain trajectories. We have shown, at least in the case of a 
separable potential, that the trajectory is required to be almost entirely along the direction of one field but that the 
orthogonal field must become more important towards the end of inflation and hence the inflationary trajectory must 
curve. In absolute terms the change to the angle of the background trajectory is small (compared to a trajectory 
which turns by a right angle during inflation), but in relative terms it must grow by at least an order of magnitude. 
This is in contrast to the previous three scenarios. 

We have reviewed two-field hybrid inflation as an explicit model of a separable potential where the conditions 
required to generate a large non-Gaussianity can be satisfied. The conditions can be satisfied for any possible com- 
bination of positive and negative 77 parameters, so the potential can be bowl shaped, a hill top or have a saddle 
point. The main conditions which must be satisfied is that the difference of the two 77 parameters must not be too 
small, rjipip — i]^^ ^ 0.1, and the value of the x field must be very subdominant to that of the ip field initially (or 
vice versa). In general, but depending on the coupling constants between the two inflaton fields and the waterfall 
field, there is a change to observables at the end of inflation, due to the fact that the surface on which the waterfall 
field is destabilised and infiation ends might not be a surface of uniform energy density. This effect is responsible 
for the inhomogcneous end of infiation scenario. It is then a model dependent question whether there will be further 
evolution to the observables during reheating in this model, this deserves further attention. 

One similarity that this hybrid infiation model has together with the quadratic curvaton scenario is that in both 
cases the initially subdominant (approximately isocurvature) field x h^^s the ratio Sx/x ~ Cx a-ppi'oximatcly constant 
and the field fluctuations do not become more non-Gaussian with time. However the effect of this light field on the 
primordial curvature perturbation grows, during inflation in the hybrid scenario and before the curvaton decay in the 
curvaton scenario. It is this non-linear transfer between the field fluctuation and C, described by the SN formalism, 
which can generate a large non-Gaussianity. Therefore the non-Gaussianity in the hybrid scenario which we have 
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studied is generated on super horizon scales during slow-roll inflation, in a similar way to which non-Gaussianity is 
generated over time in the curvaton scenario before the decay of the curvaton. The evolution of f/yr during inflation 
is explicitly calculated and plotted in [t^. This conclusion is somewhat different from that in [3J], and we plan to 
elaborate on this point in a future work. For more discussion on the distinction between non-Gaussianity generated 
by a non-Gaussian field perturbation, and non-Gaussianity generated by a non-linear transfer between a Gaussian 
field perturbation and see [93| . An example where the subdominant fields fluctuations can become non-Gaussian 
due to a large self interaction was discussed by Bernardeau (32j . 

In order to study models where slow-roll breaks down before the end of inflation, it is clearly necessary to go beyond 
a formalism based on the slow-roll approximation. We have shown how this can be done in the context a separable 
Hubble parameter instead of a separable potential and this leads to an exact expression for (the local part of) /nl in 
these models. As an explicit example an exact two-field solution with an exponential potential was given. For some 
parameter choices this leads to a strong break down of slow roll before the end of inflation, which may give rise to a 
large non-Gaussianity. Further work is also required for this model to understand how the potential may be modified 
after the end of inflation in order that reheating occurs. 

Non-Gaussianity is a topical field, in which observations have improved greatly over the last decade through both 
studies of the CMB and large scale structure. Observations so far have heavily focused on constraining the bispec- 
trum non-linearity parameter /nl- Currently the tightest constraint comes from the WMAP satellite, assuming the 
local model of non-Gaussianity this constrains the amplitude of the non-Gaussian part of the primordial curvature 
perturbation to be less than about one thousandth the amplitude of the Gaussian perturbation. This constraint is 
likely to be tightened considerably by the Planck satellite, which is currently taking data, or instead there might be 
a detection. A detection of /nl at this level would rule out simplest models of inflation, which are single field with a 
canonical kinetic term. Clearly this would be an extremely exciting result. 

However even if we are in the fortunate position of having a detection of /nl as well as improved con- 
straints/detection of the scalar-to-tensor ratio and the spectral index, there will probably still be several viable 
scenarios, as detailed in this article, which for suitable parameter choices and initial conditions can match the obser- 
vations. Fortunately non-Gaussianity is about much more than one number. The trispectrum (four-point function) 
depends on two non-linearity parameters. In general tjvl ^ (6/nl/5)^. If the current observational hints (which are 
not statistically significant) that /nl ^ 40 turn out to be true, then both the bispectrum and the trispectrum should 
be large enough for Planck to detect. Even if the bispectrum turns out to be much smaller, although for many models 
tnl is close to the lower bound, we have seen that in the model of hybrid inflation it is possible to have t^l 2> /nlj 
so the trispectrum might even be the first observational signature of non-Gaussianity. Alternatively the trispectrum 
through a large qml might give the first observational signature, as is possible in self-interacting curvaton models 
or the exact solution with an exponential potential. If /nl is detected it will also be possible to either constrain or 
detect a scale dependence of this parameter. Although it is often assumed to be constant, this is only true for certain 
simple models, and for example in the two-field hybrid inflation model it generally has a significant scale dependence. 
We have therefore seen that non-Gaussianity is an important and powerful method of constraining and distinguishing 
between the many models of inflation. 
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